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Abstract. The finite-size scaling behaviour for percolation and conduction is studied 
in two-dimensional triangular-shaped random resistor networks at the percolation 
threshold. The numerical simulations are performed using an efficient star-triangle 
algorithm. The percolation exponents, linked to the critical behaviour at corners, are 
in good agreement with the conformal results. The conductivity exponent, t' = C/ v -. 
is found to be independent of the shape of the system. Its value is very close to recent 
estimates for the surface and bulk conductivity exponents. 



1. Introduction 

Since the paper of Cardy [1] we know that, at a second-order phase transition, the 
local critical behaviour can be influenced by the shape of the system. Furthermore, in 
the case of conformally invariant two-dimensional (2D) systems, the tools of conformal 
invariance can be used, at the critical point, to relate the local critical behaviour at 
a corner to the surface critical behaviour [2,3]. The corner shape, which is scale 
invariant, leads to local exponents varying continuously with the opening angle. This 
marginal local critical behaviour has been indeed observed numerically for different 
systems (see [4] for a review) and, more recently, analytical results have been obtained 
for the Ising model [5-8]. 

In this paper, we present a numerical study of the critical behaviour of 2D random 
resistor networks with triangular shapes. This problem involves two sets of exponents, 
namely the percolation exponents and the conductivity exponents [9-11]. Although 
our main interest concerns the conduction exponents, our simulations allow us, as a by- 
product, to check the conformal predictions for the corner exponents of the percolation 
problem. 

The percolation exponents are known exactly, both in the bulk and at the surface, 
through a correspondance with the limit q — ► 1 of the g-state Potts model [12-14]. 
The conformal aspects of the critical percolation problem in finite geometries have 
been extensively studied in [15], following the work of [16]. Conformal invariance has 
been also verified in a transfer-matrix calculation of the surface percolation exponent, 
using the gap-exponent relation [17]. The critical behaviour at surface and corners 
has been considered in [18] where the conduction problem is addressed briefly. 

A recent series expansion study [19] suggests that the surface conductivity has the 
same critical behaviour as the bulk one (see [20] and references therein). The main 



1 



2 



P Lajko and L Turban 



purpose of this work is to examine, with high numerical accuracy, wether the shape of 
a finite system may have some influence on the scaling behaviour of the conductivity. 

We study the finite-size-scaling behaviour of the conductance and percolation 
probability between points located at the corners of a triangle, either on the triangular 
or on the square lattice. As in the Lobb- Frank algorithm [21], the numerical technique 
involves a succession of triangle-star and star-triangle transformations, which allows 
one to reduce the triangular resistor network to a star network in a finite number of 
steps. 

In section 2 we introduce the different correlation functions of the percolation 
and conduction problems and define the associated exponents. We also review the 
conformal results for the corner exponents. In section 3 we give a detailed explanation 
the triangle-star star-triangle algorithm. The finite-size scaling simulation results are 
presented in section 4 and discussed in section 5. 



2. Percolation and conduction correlation functions 

We consider a random resistor network for which each lattice bond has a probability 
p to have a unit conductance and 1 — p to be an insulator. Let the connectedness 
characteristic function c- • be defined as 

{1 if sites i and j are connected ,„ „ . 

(2.1) 
otherwise. 

With [• • -] av denoting a configurational average, the percolation correlation function 

Pij = Mav (2-2) 

gives the probability that sites i and j belong to the same cluster of conducting bonds. 
The average conductance given by 

Gij = \9ijU (2.3) 

where is the conductance of the system between sites i and j, plays the role of a 
correlation function, or non-local conductive susceptibility, for the conduction problem 
[22] . Let N be the number of samples taken into account in the configurational average 
and N con the corresponding number of samples for which sites i and j are connected. 
The correlation functions in equations (2.2) and (2.3) can be rewritten explicitly as 



(2.4) 



conla 



1 N N 
P„ = lim — V c% = lim -£2£ 

1 N N 1 Woo„ 

a— 1 tjUn a— 1 

Thus one can define the reduced conduction correlation function 
G 1 Ncon 

r*, = ^r= lim -T7— Y [gtfT ( 2 -5) 

13 Pij »^°°Kon 
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which gives the average conductance between two points i and j when they belong to 
the same percolation cluster. 

In an infinite system, for two points at a distance r^ = r, the correlation functions 
display a power law decay 

Pij~r- 2x r .~r-</" (2.6) 

at the percolation threshold p = p c . The exponents x and v are the scaling dimension 
of the bulk order parameter and the correlation length exponent for the percolation 
problem, respectively. The conductivity exponent £ governs the behaviour of the 
macrosopic conductivity near the percolation threshold where [23] 

£~(p-f> c )* t = vt' = (+{d-2)v P>Pc . (2.7) 

Note that C,/v = t' in two dimensions. 

As mentioned in the introduction, the percolation exponents are exactly known in 
two dimensionsthrough a correspondance with the g-state Potts model [24, 25] : 

v=\ x = ^ x s = \ (2.8) 

where x s is the scaling dimension of the surface order parameter at the ordinary 
transition. 

Recent high statistics simulations led to the following accurate estimate for the 
conductivity exponent in two dimensions [20]: 

t' = - = 0.9825 ±0.0008 C = 1-3100 ± 0.0011 . (2.9) 
v 

Low-density series expansion results [19] are consistent with a surface conductivity 
exponent £ s keeping its bulk value £. 

For a triangular-shaped system of size L, when the points i and j are located at 
corners with opening angles 9 i and • , according to finite-size scaling, one expects the 
following behaviour at criticality: 

P. = P(6-,6 ;L) ~ L- r '^' 9 ^ ~ L -<^)-<e 3 ) 

13 ' 3 (2.10) 

IV =T{6 i ,6 j ;L)~L-W^I» . 

Here x(6 i ) is the scaling dimension of the local order parameter at a corner with 
opening angle 9 { . In the following, we shall also consider the three-point correlation 
function P ijfe , which gives the probability that the points i, j and k, located at corners 
with opening angles 9 i , 9j and 9 k , belong to the same cluster. This quantity scales 
like the product of the corresponding local order parameters, i.e., as 

P ijk = P ( e i, e j,°k> L ) ^ L-iWiM - (2.11) 
at the percolation threshold. 

A dependance of the local exponents on the opening angles is generally expected 
since a wedge is a scale-invariant geometry and the angles are marginal variables for 
the local critical behaviour. The critical 2D Potts model being conformally invariant, 
one obtains the local critical behaviour for the percolation problem at a corner using 
the conformal transformation w = z e / n , which maps the half-space onto a wedge with 
opening angle 9 [2, 3]. This leads to the following expression for the scaling dimension 
of the order parameter: 

x{6) = lx s . (2.12) 




3. The star-triangle transformation 



Let us consider a finite random resistor network with the shape of an equilateral 
triangle of side L on the triangular lattice. Through a succession of triangle-star and 
star-triangle transformations, the original system can be progressively transformed 
into a star with 3 branches of length L as shown in figure 1. 

At step p, in the first part of the lattice transformation, up-pointing triangles with 
coordinates and bond conductances 9a\i,j) (a = 1,2,3) are replaced by stars 
with bond conductances 

la { l :J)— TZT " (X— 1,4,0 [0.1) 

as shown in figure 2. In the second part, down-pointing stars are transformed 
into down-pointing triangles and up-pointing triangles are relabelled as indicated in 
figure 3. Thus one obtains the bond conductances for up-pointing triangles at step 
p + 1 as 

<*-i),, ,x _ 7^(^)7^ + U) 



7i (p) (i + l,i-l)+7 2 °' ) (<,j)+73 (p) (< + l,i) 

J P+ D ( . ?) = f3 2) 

■ 2 7i (p) (i.i)+7 2 Cp) (i-l,i + l)+7 3 (p) (i,J + l) 

(P+D r 7i P) (i + l,i)7^(i,i + l) 



7< p >(i + + >y<?\i,j + 1) +7 3 (P) (< + 1,.? + 1) 
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Figure 2. In the triangle-star transformation, up-pointing triangles with conduc- 
tances g^\i,j) are transformed into stars with conductances given by 
equation (3.1). 



(i-lj+1) 



(i+lj+D 




Figure 3. The conductances g^^ 1 ^ (i, j) on up-pointing triangles at step p + 1, 
given in equation (3.2), follow from star-triangle transformations on down-pointing 
stars. The conductances of the star involved in the construction of g^ p+1 ^ (i, j) are 
indicated. They originate from different up-pointing triangles at step p. 



Note that surface bonds in up-pointing triangles at step p + 1 result from the 
transformation of incomplete down-pointing stars. The expressions given in (3.2) 
still apply in this case, provided the conductances associated with the missing bonds 
are set equal to zero, i.e., with the boundary conditions 
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Figure 4. A triangular-shaped system on the square lattice is obtained through a 
deformation of the triangular lattice when the conductances are set to zero in one of 
the principal directions. 



7 r(i,o)=o 



i = 2, L — p 



72 W (<U)=0 



j = L - p - i + 2 . 



i,j = 2,L-p 



(3.3) 



At step p = L — 1, the final star configuration is obtained after the triangle-star 
transformation has been performed. Then, for example, the conductance between 
points A and Bona sytem with size L is given by 



The same transformation, with all the conductances g$ set to zero in the initial 

configuration, can be used to reduce to a star a triangular-shaped system on the square 
lattice as shown in figure 4. 

4. Finite-size scaling results 

We have studied the finite-size scaling behaviour of the percolation and reduced 
conduction correlation functions between corners on triangular-shaped sytems. We 
worked at the percolation threshold, either on the triangular lattice (p c — 2 sin(7r/18)) 
or on the square lattice (p c = 1/2) [26]. 

On the triangular lattice (figure 1), we calculated the two-point functions 
P(^,^;L), T(^,^;L) and the three-point function P(|,|,|;L) whereas on the 
square lattice (figure 4), we studied the two-point functions P(^,j;L), P(j, f; L), 
r(f ,f ;L), T(f , f ;L) and the three-point function P(f , f , f ;L). 

The initial bond configurations were generated using two types of random number 
generators, a simple shift register algorithm and the ranlux97 generator. We checked 
that both generators led to consistent results within the statistical errors. On the 
square lattice, in order to increase the number of percolating samples, the external 
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Figure 5. Percolation on the triangular lattice: effective decay exponents for the 
two-point (a) and three-point (v) correlation functions plotted against L~ w and 
their extrapolated values (•). 




Figure 6. Percolation on the square lattice: effective decay exponents for the two- 
point (□, o) and three-point (v) correlation functions plotted against and their 
extrapolated values (•). 



bonds connecting A and B to the rest of the system in figure 4 were always assumed 
to be conducting. 

We used system sizes of the form L = 2 k up to L — 256. The star-triangle 
algorithm described in section 3 led to a computation time scaling roughly as L 2 InL. 
The number of samples generated was N = 2 x 10 7 for the triangular lattice, except 
for the largest size where N = 3 x 10 7 . On the square lattice, 4 x 10 7 samples 
were generated for all sizes. The percolation and conduction correlation functions 
P(L) and T(L) were stored as n independent averages over groups of 10 6 samples in 
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Figure 7. Conduction in triangular-shaped systems: effective decay exponents of 
the corner-to-corner conduction correlation functions on the triangular (a) and the 
square (□, o) lattices, plotted against L _ "\ The dotted lines are the best linear fits 
used for the extrapolations and the arrow indicates Grassberger's result for the the 
bulk conductivity exponent. 



Table 1. Exponents governing the decay of the percolation correlation functions 
(?7) and the conduction correlation functions (C/f) for triangular-shaped systems at 
criticality on the triangular and square lattices. The numerical values obtained for the 
percolation problem are compared to the values expected from conformal invariance. 
The last two lines give the effective correction-to-scaling exponents u> used in the 
extrapolation process. 



Triangular lattice Square lattice 



Angles (f,f) (|,f,f) (f,f) (f,f) (f.f.f) 

r) numerical 1.99 ±0.02 3.0 ± 0.3 2.01 ± 0.05 2.65 ± 0.07 3.32 ± 0.16 

T) conf. inv. 2 3 2 I T 

C/v numerical 0.9827 ± 0.0017 0.9829 ± 0.0012 0.979 ± 0.018 

uj percolation 1.05 0.76 1.14 1.06 1.13 

uj conduction 0.81 0.43 1.02 



order to evaluate the statistical errors. More precisely, as defined in equation (2.5), 
the conduction correlation function is averaged over that part of the 10 6 samples for 
which the two points are connected. 

Effective exponents at size L for the percolation and conduction problems were 
obtained using the two-points approximants 

■ hP(V2H ■nP(2L) ^. hTM-hnal ^ (41) 

The central value of the approximants are calculated by averaging P(L) and T(L) 
over all the samples. The error bars were deduced from the mean-square deviations 
a 2 of the approximants, obtained with the n statistically independent averages, as 
iyVAn-l). 
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The finite-size results for the percolation exponents are shown in figure 5 for the 
triangular lattice and figure 6 for the square lattice. The results for the conductivity 
exponents are shown in figure 7 for both lattices. 

The finite-size results were extrapolated in the following way. Assuming a single 
correction-to-scaling exponent uj, in each case we looked for the value of u> leading to 
the best linear variation at large size for the central value of the effective exponent 
as a function of L _a \ The intercept with the vertical axis gives the extrapolated 
values t] = r/ eff (oo) and £/v = £/z/| eff (oo). The statistical error on the extrapolated 
values were deduced, as above, from the mean-square deviation of the extrapolated 
exponents deduced from n statistically independent sets of approximants, taking the 
same value of u> for the linear fit. To the statistical error, we added a systematic error, 
linked to the deviation from the asymptotic regime. It was taken as the difference 
between the extrapolated values when one takes or not into account the largest size in 
the extrapolation process. Doing so, we probably overestimate the systematic error. 

Due to the non-monotonous behaviour of the effective exponent for the 
conductance between corners with opening angle j , the error bar on the extrapolated 
exponent had to be estimated differently It was obtained through an extrapolation of 
the extreme values of the effective exponents using the same method as for the central 
value. 

The extrapolated exponents for the percolation and conduction problems are 
shown in table 1. 



5. Discussion 

Let us first consider the percolation problem. The decay exponents of the two- 
and three-point correlation functions, following from conformal invariance, are easily 
obtained using equations (2.8) and (2.10)-(2.12). Our numerical results in table 1 arc 
in good agreement with the expected ones, although with a lower precision for the 
three-point exponents, due to larger statistical errors. 

The percolation probability between two corners P(8 i7 9j;L) can be written as 
the sum of two contributions. The leading one is the probability that i and j are 
connected without being connected to the third corner k which scales as L~ r, ^ i,6 '\ 
The second is the probability that the three corners are connected and it decays as 
jj-riiOifijfik) _ Thus we have a correction-to-scaling exponent equal to x(9 k ) for the 
two-point percolation probability. Another correction is due to the leading irrelevant 
operator with scaling dimension —1 at the surface [17]. This explains the value, close to 
1, of the effective correction exponent u> for two-point percolation, since the amplitude 
of the first correction is small. 

Our main result concerns the decay exponent C,/v which is equal to t' in two 
dimensions. The numerical values given in table 1 for the different geometries are all 
quite close to the value of t' for the bulk, which is given in equation (2.9). These 
results, together with the surface ones [19], strongly support the existence of a single 
conductivity scale, independent of the sample geometry The influence of the opening 
angles can be seen in the amplitudes and perhaps also in the correction-to-scaling 
exponents. 

One may notice that, apart from the case of the diagonal direction on the 
square lattice, the approximants for the conductivity exponents display a monotonous 
behaviour at large size. This is to be compared to the non-monotonous behaviour 
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observed in reference [20] for bond percolation conductivity on the square lattice. It 
makes the extrapolation of the exponents more easy, thus reducing the computational 
effort necessary for a given precision on the extrapolated values. 

Except for the diagonal direction on the square lattice, our error bars are about 
the double of the ones of Grassberger, for a maximum system size L = 256 instead of 
L = 4096 in [20]. If one accepts the universality of the conductivity exponent, which 
is strongly suggested by our results, the triangular geometry used here appears as a 
potentially efficient tool to further improve the precision on the value of Q/v. 
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